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Internal Shocks in the Magnetic Reconnection Jet in Solar Flares: 
Multiple Fast Shocks Created by the Secondary Tearing 

Instability 

S. Tanuma^ & K. Shibata^ 

ABSTRACT 

Space solar missions such as Yohkoh and RHESSI observe the hard X- and 
gamma-ray emission from energetic electrons in impulsive solar flares. Their 
energization mechanism, however, is unknown. In this paper, we suggest that 
the internal shocks are created in the reconnection jet and that they are possible 
sites of particle acceleration. We examine how magnetic reconnection creates the 
multiple shocks by performing two-dimensional resistive magnetohydrodynamic 
simulations. In this paper, we use a very small grid to resolve the diffusion 
region. As a result, we hnd that the current sheet becomes thin due to the 
tearing instability, and it collapses to a Sweet-Parker sheet. The thin sheet 
becomes unstable to the secondary tearing instability. Fast reconnection 
starts by the onset of anomalous resistivity immediately after the secondary 
tearing instability. During the bursty, time-dependent magnetic reconnection, 
the secondary tearing instability continues in the diffusion region where the 
anomalous resistivity is enhanced. As a result, many weak shocks are created 
in the reconnection jet. This situation produces turbulent reconnection. We 
suggest that multiple fast shocks are created in the jet and that the energetic 
electrons can be accelerated by these shocks. 


Subject headings: acceleration of particles - Sun: flares - Sun: corona - MHD - 
plasmas - turbulence 


1. INTRODUCTION 

In impulsive flares, strong hard X-ray emission from energetic electrons has been ob¬ 
served by Yohkoh HXT (Masuda et al. 1994; see also a review by Miller et ah 1997). The 
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origin of energetic electrons is, however, unknown. The fast-mode shock is known to be a 
possible site for particle acceleration, for example in interplanetary space (Terasawa & Sc- 
holer 1989) and supernova remnants (Blandford & Ostriker 1980; Koyama et al. 1995). To 
solve the particle acceleration problem in solar flares, Tsuneta & Naito (1998) suggested a 
diffusive shock model at the loop top. This model, however, has not yet been fully examined. 

Solar flares occur as a result of magnetic reconnection. The SXT and HXT on the Yohkoh 
satellite discovered various pieces of evidence of magnetic reconnection (Tsuneta et al. 1992; 
Masuda et al. 1994; see review by Shibata 1999). The “fast” reconnection process suggested 
by Petschek (1964) can be supported by a localized resistivity (e.g, anomalous resistivity; 
Ugai 1986; Tanuma et al. 2003; see also Shibata, Nozawa, & Matsumoto 1992; Miyagoshi 
& Yokoyama 2003). Recently, Tanuma et al. (1999, 2001) investigated the reconnection 
triggered by a shock wave from a point explosion. They found that Petschek reconnection 
occurs immediately after the “secondary tearing instability” (Furth, Killeen, & Rosenbluth 
1963; see also Magara & Shibata 1999; Shibata & Tanuma 2001), and that fast shocks are 
created in the reconnection jet. In this paper, we suggest that energetic electrons could be 
accelerated if the internal fast shocks could be created in solar flares (hgure 1). To examine 
this possibility, we study how the internal shocks are created in the reconnection jet by 
performing 2D simulations with high spatial resolution. This should remove the effect of 
numerical noise and resolve the diffusion region. 

These simulations lead us to propose that multiple fast shocks are created in solar flares 
and they are possible sites of particle acceleration. Our simulation method is described in 
the next section. In section 3, we present the simulation results. Finally, we discuss the 
results. 


2. Model of Numerical Simulations 

We solve the nonlinear, time-dependent, resistive, compressible MHD equations. A 
rectangular computation box with 2D Cartesian coordinates in the x-z plane is assumed, 
where the x and axes are in the horizontal and vertical directions, respectively. We assume 
an ideal gas, i.e., Vg = ~ l)e, where 7 is the specihc heat ratio (=5/3). The velocity and 

magnetic held are v = (n^W^) and B = Gravity is neglected. 

In the initial equilibrium conditions, a Harris current sheet is assumed, i.e., B{x,z) = 
Rotanh( 2 :/r“*)a;, Pgix,z) = Pgo + {B^/87r)[l - tanh^(^/d“*)], p{x,z) = {Wg/T) = Po + 
( 7 /To)(Ro/ 8'^)[1 — tanh^( 2 ://‘“*)], where Bq, pgQ, and po, are dimensionless variables, and 
X = (1,0). We assume in the initial conditions that the current-sheet half-thickness = 1. 
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The ratio of gas to magnetic pressure is /9 = Snpgo/BQ = 0.2 (| 2 :| 3 > The initial gas 

pressure is Pgo outside the current sheet, and Pgo + B^/Sn = Pgo{^ + 1//5) inside the current 
sheet. The total pressure is uniform. We assume Pgo = 1/7 = 0.6 and po = 1 - The sound 
velocity and temperature are Cg = {'yPg/ pY^"^ = 1 (uniform) and T = Tq = 1 (uniform). The 
initial Alfven velocity is = Bq/(A irpoY^'^ — 2-45 {\z\ 3> /'“*). We assume an anomalous 
resistivity model as follows: p = po for Vd < Vc, and p = po + a{vd/vc — lY for > Vc 
(Kliem, Karlicky, & Benz 1998; Ugai 1986; Tanuma et al. 1999, 2001), where Vd{= J/p), 
p, J, and Vc are the relative ion-electron drift velocity, mass density, current density, and 
threshold above which the anomalous resistivity sets in (Anomalous resistivity originates 
from a lower hybrid instability, ion-acoustic instability, or ion-cyclotron instability; see, e.g., 
Treumann 2001). We also assume that the resistivity does not exceed pmax = 1- In this 
paper, we assume a “background resistivity” Pq = 0.005, which is sufficiently larger than 
the “numerical resistivity” because of the grid size (Ugai 1999; Tanuma et al. 2001). The 
other parameters are a = 10.0 and Vc = 20.0. We neglect the thermal radiation and heat 
conduction. They do not affect the basic properties of magnetic reconnection (Yokoyama & 
Shibata 1997). 

We normalize the velocity, length, and time by the sound velocity (C*s), initial current- 
sheet thickness {H), and H/Cg, respectively. The units of normalization are Cg ~ 150 km 
s“^, H ~ 3000 km, and r = H/Cg ~ 20 s. The units of temperature, density, gas pressure, 
and magnetic held strength are Tq ~ 2 x 10® K, uq ~ 10® cm“^, pgQ ~ 10“^ erg cm“^, and 
Bq ~ 2 G, respectively {Bq is ~ 10 G if we assume Uq ~ 2.5 x 10^® cm“^). The grid number 
is {Ng.,NY = (13000,1300), and the grid size is uniform {Ax,Az) = (0.013,0.013). We 
assume that the top {z = -1-8.45) and bottom {z = —8.45) surfaces are free boundaries, and 
that the right {x = -1-84.5) and left {x = —84.5) ones are periodic. The simulation box size 
is {L^^LY) = (169.0,16.9). The magnetic Reynolds number is Re™* = v^^^L^/pq ~ 84500. 
We use a 2-step modihed Lax-Wendroff method. The resistivity is initially enhanced for a 
short time in the central region of the current sheet (see also Ugai 1986). 


3. Results 

Figure 2 shows the time variation of the spatial distribution of the gas pressure in the 
reconnection region. The current sheet is unstable because of the tearing instability (Furth 
et al. 1963). The magnetic dissipation time, Alfven time, and tearing instability time scale 
are ~ 200, ri“* = ~ 0.4, and ^ 9 , respectively, 

where /'“*(= 1) is the half thickness of the initial current sheet. The tearing instability is 
initiated in the current sheet by the initial perturbation. The current sheet becomes gradually 
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thinner in the nonlinear phase of the tearing instability {t ~ 7 — 18). Its length is comparable 
to the most unstable wavelength of the tearing instability, i.e., \t ~ (Magara 

& Shibata 1997; Tanuma et ah 2001), where Re™* = /tiq ~ 500. A* is ~ 25 in this 

simulation. The current-sheet thickness in this phase is k ~ (At/2)Re~y^, i.e., the current 
sheet corresponds to that of a Sweet(1958)-Parker(1957) sheet. The magnetic Reynolds 
number is Rem,t = (-^t/2)t^A‘*/j7o in this phase. Furthermore, Re^,* ~ 6000 and It ~ 0.15. 
These values explain our results well. 

The sheet becomes very long so that it becomes unstable to the tearing instability again 
at t ~ 17 — 78 (referred to as the “secondary tearing instability” in Tanuma et al. 1999, 
2001). The secondary tearing instability starts just after the Sweet-Parker sheet is created on 
a time scale of T 2 t = / VoV^a^Y^'^ ~ 0.52. The drift velocity (vd) reaches a threshold (vc) at 
t ~ 18 so that anomalous resistivity sets in. Vd, however, remains a little larger than Vc when 
18 < t < 20. Figure 3a shows that the magnetic islands are actually created at the diffusion 
region (But plasmoids are difficult to represent clearly in hgure 2a). The gaps between them 
(~ 0.3 —0.5) are not consistent with the analytic solution of the secondary tearing instability 
in this model, \ 2 t ~ 4.9Re^^2^^"^h, where Re^^^ = n™*h/?7o, ^ 2 t ~ 2.2, and Re^,^ ~ 73.5, 
respectively. The current-sheet thickness becomes l 2 t = Re^ ^'^^(A2t/2) ~ 0.047, where the 
magnetic Reynolds number is Re^ = (A2t/2)t;™*/?7o = 539. This is because the analytic 
solution is in the large Re^ approximation regime. 

The drift velocity increases a great deal above the threshold at t ~ 20 before the current 
sheet completes its collapse to a thickness of l 2 t- Because of the excitation of anomalous resis¬ 
tivity, non-steady Petschek(1964)-hke (fast) reconnection starts, and is accompanied by slow 
shocks. During fast reconnection [t > 20), the secondary tearing instability occurs in the dis¬ 
sipation region where the anomalous resistivity is enhanced. This diffusion region is well re¬ 
solved numerically. The magnetic Reynolds number is now Re^*^ = n™*/t/?7max ~ 0.38. The 
wavelength and time scale of the secondary tearing instability are A^^ ~ ~ 0.58 

and T^t = YY/ 0.0014, respectively, provided we assume that the secondary 
tearing instability starts when / ~ /j. Figure 3b also shows that the plasmoids are created at 
this time. Their wavelengths are consistent with the above theoretical value. This result can 
be explained by extension of the results of Steinolfson & van Hoven (1983), although the an¬ 
alytic solution is in the large Re^ approximation regime. Due to the bursty, time-dependent 
reconnection, many weak shocks are created in the reconnection jet. They are almost sta¬ 
tionary, when they are created outside the diffusion region, although they move inside and 
near the diffusion region. Figure ?? displays the prohles of some variables in z = 0.058, when 
t = 35.0. The prohles show that 12 weak shocks are created among 0 < a; = 6.5. Hence, 
the gaps between them (~ 6.5/12 ~ 0.54) are quite consistent with X^t (Note that the local 
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Alfven velocity is smaller than n™'* during this phase). The reconnection jet can become 
supersonic and the pressure jumps are almost standing so that they could be shocks. The 
reconnection rate and inflow velocity toward the diffusion region oscillate. As a result, the 
reconnection becomes “turbulent” (Matthaeus & Lamkin 1985, 1986; Lazarian & Vishniac 
1999; Fan, Feng, & Xiang 2004) during the late phase. In this phase, the Kelvin-Helmholtz 
instability may also occurs (see discussion). 

We examine the dependence of the results on the plasma (3 and resistivity model (? 7 o, 

Vci and a). These parameters do not affect our principal results; fast shocks are created 
in all simulations. For example, the physical values are Re™* ~ 500(/™*/1.0) (?7o/0.005)“* 
(/3/0.2)-R^ At ~ 25(/™*/1.0)^/^ (r/o/0.005)-R^ (/3/0.2)-R®, and k ~ 0.15(/™*/1.0)5/® (r/o/0.005)3A 
(/5/0.2)^R®, respectively. Furthermore, A^* ~ 0.58(/™*/1.0)^®/^^ (?7o/0.005)*®/^^ (^7max/l-0)“*^^ 
(/3/0.2)R6^ andr 2 *t ~ 0.0014(/™*/1.0)*5/*6 (7/0/0.005)**/*^ (r/max/1.0)-*/^ (/3/0.2)*R32^ if anoma¬ 
lous resistivity is assumed to be excited to ~ T^max when the current-sheet thickness is I ~ k. 

A detailed investigation of the parameter-dependence will be published elsewhere. 


4. Discussion 

We suggest that bursty, time-dependent reconnection creates many weak shocks (plas- 
moids) in the reconnection jet, which could become multiple fast shocks in actual flares. 
Reconnection produces an outflow whose velocity is comparable with the local sound speed. 
Thus, it can create fast shocks which would provide sites for Fermi acceleration. From our 
numerical simulations, the number of energetic electrons can be determined by the flux into 
the internal shocks, and is calculated to be dN/dt ~ 10^® s“* (see also Tsuneta & Naito 
1998). The reconnection site evolves in a self-similar manner (Nitta et al. 2001), so the 
evidence of bursty reconnection and internal shocks shown in this paper may be observed by 
the Solar-B satellite which will be launched in 2006. If the energetic electrons are accelerated 
in the fast shocks, the number of energetic electrons would increase with distance from the 
diffusion region. 

Our simulations show that the reconnection jet (“downflow”) produces turbulence. 
The turbulent held at the loop top helps to accelerate and confine the energetic electrons 
(Jakimiec et al. 1998). As a second important application, we suggest that the non-steady 
ejection of plasmoids can create the internal shocks. They actually appear in some of our 
simulations with different parameters, which we will examine in a future paper. This could 
explain the oscillations (Innes, McKenzie, & Wang 2003) and non-steady plasmoid ejec¬ 
tions (McKenzie & Hudson 1999; Asai et al. 2004) in the downflows. In fact, some results 
with different parameters show the oscillation and multiple oblique shocks in the jet, which 
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could be due to the Kelvin-Helmholtz instability (Biskamp, Schwarz, & Zeiler 1998; Tanuma 
2000; Arzner & Scholer 2001; Tanuma & Shibata 2002, 2003b, 2004). Furthermore, spatial 
patterns that are similar to the multiple oblique shocks appear in galactic jets (e.g., Hardee 
et ah 1991). 

The non-steady ejection and bow shock propagation can also play an important role in 
slowly drifting structures (Karlicky 2003), and narrowband dm-spikes (Barta & Karlicky 
2001) in the upward jet. The time scale for the tearing instability in the diffusion region 
is ~ 30 (//lO^ cm)^/^ (ua/ 10® cm s“^)“^/^ (T/10® K)^/"^ s if we assume Spitzer conductivity 
and that the diffusion region thickness is 1. This explains the time scales associated with the 
above phenomena (Karlicky 2003). Our results can also be applied to particle acceleration 
in the Galactic plane (Tanuma & Shibata 2003a) etc. 
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Fig. 1.— Schematic illustration of the internal shocks in the reconnection jet in solar flares. 



Fig. 2.— The gas pressure. The drift velocity reaches the threshold at t ~ 18. Petschek-like 
reconnection starts at t ~ 20. Bursty, time-dependent fast reconnection creates many weak 
shocks in the reconnection jet. 



Fig. 3.— (a) Bz and B vectors in diffusion regions. The white circles indicate the large 
Bz regions. It shows that small magnetic islands are created, although their gas is smaller 
than the theoretical wavelength of the secondary tearing instability {X 2 t ~ 2.2). (b) The 
gas pressure and v vectors at the diffusion region. The gaps between the shocks are quite 
consistent with the theoretical value ~ 0.58; shown by the arrow). 
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Fig. 4.— The profiles of some variables in 2 : = 0.032 at t = 35.0 between ’a’ and ’/?’ plotted in fig 
2 . Gray boxes show the pressure jumps created in the reconnection jet. They are almost standing 
shocks. The Rankine-Hugoniot relation is satisfied, for example, by more than 91 % at the jump 
named ’B’. In the 3rd plot, the local sound velocity (Cg) is also displayed by dashed-dotted lines. 
The velocity of jet is comparable with ICsl so that the jet can be supersonic and can be weak 
shocks. The pressure jumps are resolved by more than 10 grids. The bottom plot is same with 3rd 
one but in 2 : = 0.025 at t = 60. labelprofile 
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